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Abstract Since its discovery as an X-ray source with the Einstein Observatory, 
the hot X-ray emitting interstellar medium of early-type galaxies has been studied 
intensively, taking advantage of observations of improving quality performed by 
the subsequent X-ray satellites ROSAT, ASCA, Chandra and XMM — Newton, and 
comparing the observational results with extensive modeling by means of numeri- 
cal simulations. The hot medium originates from the ejecta produced by the normal 
stellar evolution, and during the galaxy lifetime it can be accumulated or expelled 
from the galaxy potential well. The main features of the hot gas evolution are out- 
lined here, focussing on the mass and energy input rates, the relationship between 
the hot gas flow and the main properties characterizing its host galaxy, the flow be- 
havior on the nuclear and global galactic scales, and the sensitivity of the flow to 
major galaxy properties as the shape of the mass distribution and the mean rotation 
velocity of the stars. 



1 Introduction 

X-ray observations performed during the late 1970's first revealed that early-type 
galaxies (ETGs) emit soft thermal X-rays and then host an interstellar medium 
(ISM), of a mass up to ~ 10 10 M Q , that had not been discovered previously at other 
frequencies because of its high temperature (T ~ 10 7 K; see Fabbiano, this volume). 
Such a medium was the long sought phase fed by stellar mass losses and predicted 



by the stellar evolution studies (e.g., 11421 17611 ): in fact, both isolated ETGs and those 



in groups or clusters emit soft thermal X-rays ll40ll . which provides simple evidence 
that the hot ISM is mostly indigenous rather than accreted intergalactic medium. Be- 
sides solving a puzzle, this discovery opened the study of an important component 
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of ETGs: it is a major ingredient of galactic evolution, see for example its role in 
feeding a central supermassive black hole and maintaining an activity cycle and star- 
formation (e.g., 144147116511 . and Ciotti & Ostriker, this volume), or tha t in polluting 
the space surrounding ETGs with metals, via galactic outflows ( jl06ll . and Pipino, 
this volume), or finally in responding to environmental eff ects, as interaction with 



neighbors, stripping, sloshing, and conduction (e.g., 11L I61LI115I1 . and Sarazin, this 
volume). 

The most striking X-ray observational feature of ETGs is the wide variation in 
their luminosity (Lx) values, of >2 orders of ma gnit ude at any fixed galactic optical 
luminosity Lb, when Lg > 3 x 10 10 Lb © 18711 . This feature cannot be ex- 



plained by distance uncertainties, since a variation of the same size is present even 
in the distance-independent diagram of Lx/Lq versus the central stellar velocity dis- 
persion a c [e.g., [39(]. A number of reasons have been proposed as responsible for 
this large variation, of environmental nature (see Sarazin, this volume) or linked to 
the possibility for the gas content to evolve substantially during the ETG lifetime. In 
this latter context, many studies investigated with numeric al simulat ions the dynam- 
ical evolution of the hot ISM in ETGs IHSIIlSliSEMIIIi- In more recent 



times, the effect of feedback from a central supermassive black hole (MBH) has 
revealed as another potential contributor to the variation of the hot ISM luminosity 
(see Statler, and Ciotti & Ostriker, this volume). 

Below I briefly review our current knowledge about the feeding and the ener- 
getics of the hot gas flows, concentrating on their dynamical state as a function of 
galactic mass and other major galaxy properties. Section |2] gives an updated sum- 
mary of the fundamental elements entering the problem, the mass and energy inputs 
to the flow; Sect.[3]presents the general case for the evolution of the flow, pointing 
out the different behavior on the global and nuclear galactic scales, and making also 
use of a representative ETG; the effect of a central MBH, acting as a gravitating 
point mass, and the expectations for accretion feedback, are presented in Sect. HJ 
finally, Sect. |5]discusses the observed sensitivity of the flow to major galaxy prop- 
erties as the flattening of the mass distribution, the mean rotational velocity of the 
stars, and the shape of the stellar profile. 



2 Feeding and Energetics of the Hot Gas Flows 

In this Section the fundamental processes and quantities at the basis of the origin 
and evolution of hot gas flows are introduced: their feeding via stellar mass losses 
(Sect. 12. It , their heating via type la supernovae explosions (Sect. 12. 2t , and their 
energy budget (Sect. |2.3b . 
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2.1 The Stellar Mass Loss Rate 

In ETGs the gas is lost by evolved stars mainly during the red giant, asymptotic gi- 
ant branch, and planetary nebula phases. These losses originate ejecta that initially 
have the velocity of the parent star, then interact with the mass lost from other stars 
or with the hot ISM, and mix with it. The details of the interaction are controlled 
by several parameters, like the velocity of the mass loss relative to the hot phase, or 



the density of the ambient ISM 17711 . Parriott & Bregman ([89] and [12]), modeling 
the interaction with two-dimensional hydrodynamical simulations, found that most 
of the continuous mass loss from giant stars is heated to approximately the tem- 
perature of the hot ISM within few parsecs of the star; in the case of mass ejected 
by planetary nebulae, about half of the ejecta separates and becomes hot, and the 
other half creates a narrow wake that remains mostly cool, unless turbulent mixing 
allows for its heating on larger scales. Far infrared observations allow us to measure 
directly the stellar mass loss rate for the whole galaxy (M*); this was for example 
derived for nine local ETGs from ISO data yfl. When rescaled by the luminosities 
of the respective galaxies, the values of M* were found to vary by a factor of ~ 10, 
which was attributed to different ages and metallicities. The average of the observed 
rates was M* = 7.8 x 10~ 12 Lb{Lb,q) M yr _1 , and was found to be in reasonable 
agreement with previous theoretical predictions j^]. 

According to single burst stellar population synthesis models, most of the stellar 
mass is lost at early times, before an age of ~ 2 Gyr. For example, the mass lost by 
stars at an age of 2 Gyr is 25% of the total initial mass, for the Salpeter initial mass 
function (IMF), and 36% for the Kroupa IMF; at an age of 12 Gyr, the additional 
loss is of ~ 4% and 6% of the initial mass, for the same two IMFs respectively J75I1. 
Figure Q] shows the trend of M*(f) with time, estimated from the models in 117511 for 
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Fig. 1 Evolution of the stellar mass loss rate M,(f) according to the models in |75|| . for a Salpeter 
IMF (solid line with best fit of Eq.[T]in red) and for a Kroupa IMF (dotted line with best fit of Eq.fT] 
in blue). Solar abundance is assumed, and a total stellar mass of 1 Mr. at an age of 12 Gyrs. 
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solar metal abundance. At an age of >2 Gyr, this trend can be approximated as: 

M,(t) = IQ- 12 A xM*(M ) ff 2 L3 (Moyr- 1 ), (1) 

where M* is the galactic stellar mas^] at an age of 12 Gyr, tu is the age in units of 
12 Gyrs, and A = 2.0 or 3.3 for a Salpeter or Kroupa IMF (see Fig.Q}. 

The relation above agrees well with previous theoretical estimates I201 17(1 . Tak- 
ing the stellar mass-to-light ratio in the B-band at an age of 12 Gyr (M*/Lg =9.17 
and 5.81, respectively for the Salpeter and Kroupa IMFs, [75]), Eq. Q] gives M*(12 
Gyr)= B x 1 0" 1 1 L B (Lb,® ) M Q yr~ 1 , with B= 1 . 8 or B= 1 .9 for the Salpeter or Kroupa 
IMF. The latter relation gives a rate that is roughly double as large as the average 
of the observational estimates quoted above (f 3]), that however has a large varia- 
tion around it, partly explained by differences in the ages and metallicities of the 
observed ETGs. 



2.2 The Type la Supernovae Mass and Energy Input 

The total mass loss rate of a stellar population M is given by the sum M(t) = 
Mt,(t) +Msn(0' where one adds to M*(f) (discussed in the previous Sect. 12. It the 
rateMsN(f) of mass lost by type la supernovae (SNIa) events, the only ones observed 
in an old stellar population [e.g., [111]. The mass input due to SNIa's is Msn(0 = 1.4 
M© 7? sn (f ) M© yr~ 1 , where Z?$n (t) (in yr~ 1 ) describes the evolution of the explosion 
rate with time, since each SNIa ejects 1.4 M©. In a detailed scenario for the SNIa 



precursors and their subsequent explosion past a burst of star formation 115 ll 15211 . 
^sn(0 experiences a raising epoch during the first 0.5-1 Gyr, at the end of which it 
reaches a peak, and then decreases slowly with a timescale of the order of 10 Gyr, 
and accounts for the present day observed rate. A parameterization of the rate after 
the peak, in number of events per year, is 

«sn(0 = 0.16(// /70) 2 x 10- 12 L B (L B , e ) ff/ (yr" 1 ), (2) 

where Hq is the Hubble constant in units of km s _1 Mpc~\ Lb is the present epoch 
galaxy luminosity, and s describes the past evolution; when t\2 = 1, Eq. [2] gives 



the rate for local ETGs in H18H . that has an uncertainty of ±30%. Recently, new 
measurements of the observed rates of supernovae in the local Universe, determined 
from the Lick Observatory Supernova Search (LOSS; lr69lo . gave a SNIa's rate in 



ETGs consistent with that in 1 1 811 . For the rate in Eq. |2] and for Hq = 70 km s~' 
Mpc~', one obtains Msn(12 Gyr) =2.2x 10~ 13 L B (L BiQ ) M yr~', that is almost 
~ 100 times smaller than the "normal" stellar mass loss rate M„ (12 Gyr)w 2 x 
10~ n L b (L b ,q) M yr~ ! derived above (Sect.Ej). 



1 The stellar mass M t changes very little at late epochs, for example by < 1% for a variation of ±2 
Gyr at an age of 12 Gyr. 
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SNIa's provide also heating (Sect. [2~3l below) at a rate Lsn(0 that is the product 
of the kinetic energy injected by one event (£sn ~ 10 51 erg) times the rate /?sn(0' 
This assumes that the total of £sn is turned into heat of the hot ISM, an assumption 
that is clearly an overestimate, but not totally unreasonable for the hot diluted gas 
(S). Then 

LsnW =Esn xR SN (t) =5.1(H /70) 2 x lO 3O L fi (L B , ) ff 2 5 (ergs" 1 ), (3) 

and is plotted in Fig.[2]for ?i2 = 1 and Hq = 70 km s Mpc~'. The SNIa's specific 
heating, given by the total SNIa's heating per total injected mass (approximated 
hereafter with M t ) is 



|K = 1 . 6( « o/7 0)W»tf-.(er S M-) 



where M* given in Eq. Q]has been used. At an age t\2 = 1, for Hq = 70 km s 1 
Mpc" 1 , and using M*( 12 Gyr) ps 2 x 10~ n L B (L B ,e) M yr -1 (valid for both IMFs, 
see below Eq.[TJ, one gets Lsn/M* ~ 8 x 10 48 erg Mq , that is a significant heating 
(with respect to, e.g., the specific binding energy of the gas, see Sect. l2.3l and Fig.|2]i. 

During the galaxy lifetime, Lsn/M* x fl 3 ~ s < an d then can increase or de- 
crease with time increasing, with consequences on the secular gas flow behavior 
(1251 28 , 70ll : Sect. [3}. Early hydrodynamical simulations of hot gas flows enlighted 



the importance of the cosmological evolution of the SNIa's rate to avoid exces- 
sive, and unobserved, mass accumulation at the galactic centers: if L$n decreases 
faster than (s > 1.3), then at early times it can be large enough to drive the 
gas lost by stars in a supersonic wind, that later can become subsonic and evolve 
into an inflow ( ll20Tl : Sect. l2.3l below). Subsequently, mainly following Chandra and 
XMM —Newton observations, it was realized that another major source of ISM heat- 
ing can be provided by the central MBH. The MBH heating is not sufficient by itself 
to avoid long-lasting and massive inflows at early times, but when coupled with the 
SNIa's heating, it can help the galaxy degassing and prevent large mass accumula- 
tion at the galactic centers, independently of the relative rates of M. t (t) and /?sn(0 
(see Ciotti & Ostriker, this volume). Recent estimates of the slope s agree with a 



value around s - 1 ( El III III El). 



2.3 Energetics of the Gas Flows 

The material lost by stars is ejected at a velocity of few tens of km s _1 and at a 
temperature of <10 4 K ( 18911 ): it is then heated to X-ray emitting temperatures by 
the thermalization of the stellar velocity dispersion (as it collides with mass lost 
from other stars or with the ambient hot gas, and is shocked) and of the kinetic 
energy of SNIa's events (Sect. 12. 2t . The first process heats the ISM at a rate 
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Lrj = — 



1 M*(0 



2 M* 



47tr 2 p t (r)o 2 (r)dr, 



(5) 



where p*(r) is the stellar density profile, and a(r) is the trace of the local stellar 
velocity dispersion tensor. The latter can be obtained by solving the Jeans equations 
for an adopted galaxy mass model (e.g., 17[]), and for an assumed stellar orbital 
anisotropy. So doing, one derives that the stellar heating L a is a few times lower than 
that provided by SNIa's, L$n, for reasonable stellar and dark matter distributions 
( l20ll : see also Sect. EDI 

In case of gas flowing to the galactic center (inflow), as the gas falls into the 
potential well, power is generated that can heat the gas; for a steady inflow of 
through the galactic potential down to the galactic center, this power is given by 



■-,1 

0J 10" 



-■J- 
o 



/+ = 

^grav 



(6) 







• 

• • 


1 r 














» ^ - ^. - ' 

^ ^ " ^ *-* 
^ " ,,. * ^ " 


L - ' 
- B XRB| 

^ ^ ^ ^ - 


I v 




^ - " * " ; 





log r. B (t, bg ) 



Fig. 2 The X-ray luminosity of the hot ISM Lxjsm from Chandra observations, versu s the galactic 
Lb (from corrected apparent magnitudes from Hyperleda, and distances from J34LI 1 20tl . for Hq = 70 
km s -1 Mpc~'). Lxjsm is from l34h (circles), l30ll (triangles), |82|| (squares); downward (upward) 
aiTows indicate upper (lower) limits. The dashed lines show the mean total luminosity from X- 
ray binaries derived from Chandra data (in blue), and its ± 1 o~ uncertainty, from floll . All X-ray 
luminosities have been recalculated for the 0.3-5 keV band, assuming the spectral shape in the 
original references. The solid line is L^jv (Eq. [5}; the green line is LJ rav (Eq. |7J, calculated for 
M„(12 Gyrs) of the Kroupa IMF (Sect. [2Tt , and representative galactic mass models (see also 
Sect. [5J with: a stellar de Vaucouleurs profile, a c from the Faber-Jackson relation, an effective 
radius R e from the Fundamental Plane relation ((J), a NFW dar k mass pr ofile with a dark-to- 
luminous mass ratio of 5 in total (|63ll), and of 0.6 within R e ([TJ, [Til El), a ratio of the scale 
radii of the dark and stellar mass distributions of = 1 .5. 
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where (r) is the total potential [and the equation above applies to mass distribu- 
tions with a finite value of (0)] . In outflows, the work done against the gravitational 
field to extract steadily and bring to infinity the gas shed per unit time M*(f) is 



M, 



L zra V = ~ '"^ /. 47rr 2 p* (r ) {r)dr. (7) 







Lg rav is the minimum power required to steadily remove the stellar mass loss, be- 
cause of energy losses due to cooling, but these are expected to be low in the low 
density of an outflow. Figure [2] shows an example of L7 rav as a function of Lb for 
representative two-component galaxy mass models. 

Even though stationary conditions are unlikely to be verified, the quantities above 
are useful to evaluate in first approximation the energy budget of the flow, and then 
predict its dynamical state and the gas content of an ETG. In a more compact nota- 
tion, the integral in Eq.[7]can be expressed as 

L- rav =M*(t)G?r-(^,p) (8) 

where a c is the central stellar velocity dispersion, and the dimensionless func- 
tion r depends on the depth and shape of the potential well, via the variables 
Sf,=MhjM % (M/, being the total dark halo mass), and j3 = ru I r„ (r% and being the 
scale radii of the dark and stellar mass distributions; see j20ll2lll92ll for examples of 
r for vari ous s tellar density profiles), F increases for increasing and decreas- 
ing j3 (e.g., llOOll f 2 ! For reasonable galaxy structures, F does not largely vary with 
and j3, as for example empirically demo nstrated by the existence of scaling laws 



as the Fundamental Plane of ETGs ( 110511 ). Equation [8] then shows how the larger 



a c , the harder for the gas to leave the galaxy; for example, per unit gas mass, one 
has L^^lM* °= a c 2 . L~ av is expected to steeply increase with Lb, since on average 
the larger is o c , the brighter is the optical luminosity of the ETG (from the Faber- 
Jackson relation Lb « <7 C 4 ), and since also M* is proportional to Lb (Sect. 12. lb . For 
a fixed galaxy structure [i.e., a fixed F~(^,/3)], and using the Faber-Jackson rela- 
tion, one has L~ v °= L l B 5 , a trend close to the green line in Fig.|2](that has a slope 
of ~ 1 .8, though, since it derives from the full Lb — O c relation for ETGs, that has a 



shallower slope at the lower Lb, Bill ). Since the available heating to extract the gas 
(L$tv) increases with Lb as well, it is useful to evaluate the run of the ratio between 
the power required to extract the gas (L~ v ) and that given by SNIa'^3 



L g rav (t) s _ l: 



oc^- 13 c7 c 2 r-(^,j3). (9) 



2 An expression equal to that in Eq.[8]can be written for L+ v , just replacing F with a function 
r + (3i,P)\ the latter has the same trend as F ~ to increase for increasing iff, and decreasing j8; for 
reasonable galaxy mass models, F + ~ 2F~ ( llOfJl l. 

3 In order to account for all the heating sources, one should add to the denominator of Eq.|9]also 
L(j, that has been neglected for simplicity, since typically L « L$n, as written below Eq.[5] 
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In a first approximation, (most of) the galaxy will host an outflow if this ratio has 
always been lower than unity, and an inflow soon after it becomes larger than unity 
( I20I0 . The time evolution of the ratio is determined by the value of s — 1 .3; recent 
progress indicates s>l (end of Sect. 12. 2b . which produces a ratio in Eq.|9]decreasing 
with time. This means that with time increasing the gas has a tendency to become 
hotter and, if outflowing regions are present, the degassing becomes faster (see also 
Sect. 13.3b . Equation |9] also indicates the underlying cause of the average Lx.ism — 
Lb correlation (Fig. that is the increase of L~ rav /LsN with , and then with 
Lb (provided that the function T does not vary widely with Lb, as expected for 
reasonable galaxy mass models). 

The relative size of LJ rav and Lsn can be estimated from Fig. [2] L~ mv < L$n for 
Lb^3 x 10 10 Lb.o, therefore in these ETGs we can expect outflows to be important, 
and then low Lx values. This "prediction" has been confirmed to be true recently, 
thanks to Chandra observations of low Lb galaxies, probing for the first ti me g as 



emission levels even below those of the X-ray binaries emission 

(E3 M H3 SSI]). 



For Lb > 3 x 10 Lb,©, instead, L~ v > Lsn an d the SNIa's heating is insufficient to 
prevent (at least some) inflow. At high Lb, however, a very large variation of Lx is 
observed, from values typical of winds to values even larger than predicted by mod- 
els for global inflows in isolated ETGs (see Sect. |3]below). This has a few possible 
explanations: on one hand, there is the high sensitivity of the gas behavior to varia- 
tions in the parameters entering Eq.|9](the dark and stellar mass, their distribution, 
the orbital structure, and possibly also Lsn, can all vary at fixed Lb), as discussed in 
Sects. [3] and [5] below (see also Il20l l92ll ). On the other hand, the simple arguments 
above do not consider important factors that can influence the hot gas content, as in- 
jection of energy from the nucleus (see Ciotti & Ostriker, this volu me: I198|| ) and/or 
environmental effects (Sarazin, this volume; lfl4l [l5 , 55 , 58 , 93 , 1 15 , ll KjI p. The 



considerations in this Section can account for the average trend of Lx.ism with Lb, 
but significant effects can be superimposed by these factors. 



3 Decoupled Flows and variation in 



ISM 



See the full chapter (chapter 2 in Hot Interstellar Matter in Elliptical Galaxies, 
Springer, 2012; http://www.springer.com/astronomy/book/978-l-4614-0579-5) 



3.1 Gas Temperature and Galaxy Structure 

In the simulations described above (Sect. |3), the average emission weighted gas 
temperature T gas ranges between 0.3 and 0.8 keV, for a large set of galaxy models 
with different Lb, dark matter fraction and distribution, and SNIa's rate. This range 
of T gas compares well with that of the gas temperatures recently determined using 
Chandra data (e.g., lflo[l3~5 , 82]). Both in the observational results and in the models 
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T gas shows a trend to increase with o c . For example, the statistical analysis of the 
gas temperature for a sample of luminous ETGs observed with ROSAT indicates a 



correlation of type a c « 7™ 6±a , although with a large degree of scatter about this 



fit ( 08811 ). Both the size of the T gas values and their trend with <J C behave as expected, 
since the gas temperature cannot be much different from the virial temperature T V ( r 
of the galaxy potential well. In fact T v i r is defined as 

T vir =±^ [4nr*p*(r)o 2 (r)dr (10) 
3k M* J 

where k is the Boltzmann constant, jim p is the mean particle mass, m p is the proton 
mass, and a(r) is the three-dimensional velocity dispersion (as in Eq.[5]l. As already 
done for Lj~ rav and L grav in Sect. 12.31 and in analogy with Eq.|8] r vi , can be expressed 
as T V i r = nm p d£Q{ffl,fS)/k, with Q. < 1. r vi > is then proportional to a]:, which 
explains the trend of T gas with a c present in the models, and is close to the trend 
shown by the observations (HH). A simplified version of the virial temperature in 
Eq.[lO]is often used, i.e., T a = \im p a}. /k; this somewhat overestimates the true T v ,>, 
since Q < 1 . One can notice that T v f r is also the temperature linked to the gas heating 
provided by the thermalization of the stellar random motions (Eq. [5}. Therefore the 
values of r v; > are expected to represent a lower boundary to the values of T gas , due 
to the importance of additional heating mechanisms (as that provided by SNIa's). 
(abridged) 



3.2 Reasons for Decoupling 



See the full chapter (chapter 2 in Hot Interstellar Matter in Elliptical Galaxies, 
Springer, 2012; http://www.springer.com/astronomy/book/978-l-4614-0579-5) 



3.3 The Gas Flow in a Tes tease ETG 



This Section presents the evolution and the properties of the flow for a testcase 
ETG, whose optical properties place it where the Lxjsm variation is of ~ 2 orders 
of magnitude (Fig. 0: Lb = 5 x 10 10 Lb and <J C = 260 km s _1 , from the Faber- 
Jackson relation; the stellar mass profile follows a Sersic law with index n = 5, as 
appropriate for the chosen Lb (e.g., l65ll ): the effective radius R e = 6.5 kpc, from 
the Fundamental Plane relation (15[]). The dark halo has a NFW profile with &=4 
and a scale radius larger than that of the stars (j3 = ri,/R e = 1 .5). 
(abridged) 

The gas velocity, temperature and density profiles at representative epochs are 
shown in Fig.[5J [...] the 0.3-8 keV surface brightness profile at the present epoch 
is shown in Fig. 0] As typical of steep mass models, a central inflowing region 
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is present from the beginning; due to the secular increase of the specific heating 
of the gas (Eq. |4] Sect. 13.2b , with time increasing the infall velocity decreases in 
modulus, the stagnation radius migrates inward (Fig. |3). By the present epoch a 
quasi-stationary configuration establishes, by which the gas leaves the galaxy at a 
rate almost equal to the rate M* at which gas is injected by stars, and the small 
difference goes into the central sink. 

The final gas temperature profil e [. . .1 d ecreasing outward is common among 
ETGs observed with Chandra (e.g.. |35ll49[ 17810 : other profiles often occurring are 
flat or outwardly increasing ones ( 13511 . see also Statler, this volume). To obtain the 
latter two shapes, other ingredients are required with respect to those included here, 
as AGN feedback depositing heating outside the central galaxy core, thr ough a j et 
or rising bubbles (e.g., SEE 53 El]), or as an external medium (SH-ImEMI). 

(abridged) 



3.3.1 Variations in the Testcase ETG and in L X jsm 



See the full chapter (chapter 2 in Hot Interstellar Matter in Elliptical Galaxies, 
Springer, 2012; http://www.springer.com/astronomy/book/978-l-4614-0579-5) 



Hot gas flows 



11 



10- 13 


■- 1 


i i ' ' ' 4 


10-" 










r\ Q Q 1,_17 

u.j— o Kev 


io-"> 


r \ 




10-M 


— ■ ^ 


. ^ \ = 


10-" 




^\ \ X : 
\\ \ \ 1 


io- |B 




\\ V. : 




\\ \\ i 


10-'» 




\\\\ : 




\\\ 1 


10-" 




\\\ \- 






\S\ ■ 


io- el 


r i 


\\V 



0.1 1 I"' 

R(kpc) 

Fig. 4 The surface brightness profile for the model in Fig.[3](blue line), and for two other models 
whose hydrodynamical quantities are shown in the next Fig. [5] with the same line type, all at an age 
of 13 Gyr. The green line follows the optical profile of the testcase ETG, and is normalized to the 
emission of unresolved X-ray binaries, as 20% of their total collective luminosity (the unresolved 
fraction is < 25% for a sample of local ETGs observed with Chandra, excluding the very hot gas 
rich ones; floll V 



4 The Nuclear Scale 



Accretion to the center is commonly present, though from a small region, for the 
models in Fig.[5j therefore, this Section explores what modifications to the flow are 
expected from the addition of a central MBH. [...] It is now clear that MBH feedback 
is unavoidable on cosmological timescales; however, for timescales much shorter 
than the cosmological one, and closer to the present epoch, it is not fully understood 
yet how it works in detail. It is then interesting to consider a few basics aspects of 
the flow behavior for models with a central MBH but neglecting feedback, as: how 
the flow is affected by the MBH gravity, the similarity of the accretion flow with a 
Bondi flow, how much accretion energy is expected, and how the picture outlined in 
the previous Sections is modified. 



4.1 Gravitational Heating from the MBH and Comparison with the 
Bondi Accretion 

The addition to the mass model of the testcase ETG (Sect. [33]) of a central gravitat- 
ing mass Mbh = 3.6 x 10 Mm, as predicted by the Magorrian relation, that remains 
constant with time, produces a flow that is not influenced by the MBH in the outer 
region; within the central ~ 100 pc, instead, the infall velocity and the temperature 
become larger (Fig. [6]), the density lower, and the mass inflow rate to the center M;„ 
smaller (from 0.12 to 0.06 M©yr _1 at a radius of 5 pc). The most relevant differ- 
ence caused by the addition of the MBH is the higher gas temperature at the center; 
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R (pc) 

Fig. 5 Radial trend of the main hydrodynamical quantities at a galaxy age of 13 Gyr, for the 
testcase ETG in Fig.[3j(blue lines), and for three galaxy models with the same Lg and cr t . except 
for the following differences: a SNIa's rate increased by 20% (dashed line); n = 6, M t /Lg = 6.1 
(solid line); n = 4, M t /Lg = 8.5 (dotted line). The corresponding gas luminosities are respectively 
Lxjsm = 1.1 x 10 39 , 4 x 10 38 , 3.0 x 10 40 erg s _1 . The brightness profiles of the models shown with 
a dashed and solid line are plotted in Fig.|4](see Sect. [53}. 



this is due to the heating of the gas produced by the sharp increase of the stellar 
velocity dispersion within the sphere of influence of the MBH (of radius of the or- 
der of GMbh/<3c = 20 pc, where a c is the velocity dispersion without the MBH; 
Fig.|7]i, and by the compression of the gas caused by the gravity of the MBH. The 
MBH heating succeeds in "filling" the central "hole" in temperature of the model 
without MBH, and it may even be able to create a small central peak in the gas 
temperature (Fig. |6). The difference in central temperature takes more time to es- 
tablish for models with larger central gas density (as that with n = 4 in Fig. 13), 
and it may even not take place within the present epoch. In models where the flow, 
without the MBH, keeps hot down to small radii (as that of the black solid lines in 
Fig-0, the change in temperature is far less dramatic (an increase by 25% of the 
central temperature), and M,„ remains substantially unchanged at very small values 
(M in ~ 1(T 4 - lO^Msyr" 1 ). 

Another important property of the flow is the value of its velocity with respect 
to the sound velocity v s , calculated for example for the adiabatic case (y = 5/3; 
v v = J ykT I jJ.m p ), and also shown in Fig. [6] Without a MBH, in cases of "cold" 
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Fig. 6 Left: the velocity and temperature profiles for the testcase ETG of Sect. I3.3l (in blue), and 
for the same galaxy model with a central MBH of mass 3.6 x 10 8 M; ; . added (in black), at an age of 
13 Gyr; the dashed lines show — v s for y = 5/3, for the two models respectively. Right: the cooling 
time and its ratio with the infall time, for the models on the left. See Sect.|4]for more details. 



accretion as for the testcase ETG, the flow becomes supersonic close to the galac- 
tic center; with a MBH, accretion is hot and keeps subsonic, while the flow tends 
to reach the sound velocity at the innermost gridpoint (that, for this series of runs, 
has been put to 1 pc). In fact, in absence of momentum feedback, it is unavoidable 
for the flow to tend to the free fall velocity close to the MBH, where the potential 
energy per particle becomes larger than the thermal energy (and this is reproduced 
by the inner boundary condition of a vanishing thermodynamical pressure gradi- 
ent). The MBH heating also causes the cooling time t coo i to become much larger 
than the inflow time f;„y within the central ~ 100 pc (Fig.|6]l. Both properties (the 
inflow velocity that tends to v s , and t coo [ >> f in /) characterize also the Bondi (1952) 
solution for spherically symmetric accretion on a central point mass, from a nonro- 
tating polytropic gas with given density and temperature at infinity, in the adiabatic 
case (y = 5/3). This fact provides some support to a com monly used procedure to 



estimate the MBH mass accretion rate of ETGs (e.g., 11711,11 1411 ). that is the use of 



the analytic Bondi (1952) formula, replacing infinity with a fiducial accretion radius 
race = 2GMbh/v% ( 04811 ). and calculating v s as close as possible to the MBH. This is 
not a trivial aspect since there are additional ingredients in the galactic flow that are 
not included in the Bondi (1952) analysis, but are accounted for by the simulations, 
as: 1) the presence of mass and energy sources, as the stellar mass losses and the 
SNIa's heating; 2) the possibility of cooling; 3) the fact that r acc is not a true in- 
finity point, since the gas experiences a pressure gradient there; 4) the contribution 
of the galactic potential added to that of the MBH. The simulations however have 
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Fig. 7 The projected stellar velocity dispersion profile for the galaxy model of the testcase ETC 
(Sect. [331 dotted line), and for the same model with a central MBH of mass 3 x 10 8 M. : discussed 
in Sect.|4](solid line). Both models have an aperture velocity dispersion of 260 km s within R t ,/8. 



some limits too: for example, the discrete nature of the stellar distribution becomes 
important where the accretion time on the MBH (f,-„y ~ 10 4 — 10 s yrs from 10 pc, 
in the simulations) is comparable to (or lower than) the time required for the stellar 
mass losses to mix with the bulk flow H3l), or to the time elapsing between 
one SNIa event and the next (see also 111 1911 ). Another limit is that some form of 
accretion feedback is also expected, as briefly outlined in the next Section. 



4.2 The Importance of the Energy Output from Accretion 

The mass inflow rates of the models in Fig. [5] with a central MBH added, range 
fromM,„ = 4 x 10~ 4 (model with solid line) to 4 x 10 _1 M yr _1 (model with dot- 
ted line), at an inner gridpoint of 5 pc, at the present epoch. If totally accreted, the 
largest of these M,„ releases an accretion power L acc ~ 0.1M,„c 2 = 2 x 10 45 erg 



s ( 14811 ). a large value that can have a significant impact on the surrounding hot 
ISM, depending on the fraction that can interact with the ISM and be transferred 
to it. Lacc is mostly in radiative form at high mass accretion rates; more precisely, 
the radiative efficiency of the accreting material can be written in a general way 
as e = 0.2 x 100m/(l + 100m), where m = Mi n /M£dd is the Eddington-scaled ac- 



cretion rate, mdM Edd =L Edd /0.lc 2 = 2.2 x 1O~ 8 M B// (M ) M yr _1 (|25|]). When 
m>0.01, then e - 0.1 - 0.2, and L acc is mostly in radiative fo rm; when m <SC 0.01, 
then e ~ 20m, as for radiatively inefficient accretion flows ( l84ll ). and the radia- 
tive output becomes negligible. In the low-e regime, the output of accretion may be 
dominated by a kinetic form (J3, S 62. 13]). For example, from the energy input 



by accretion feedback to the hot coronae of a few nearby ETGs, with M,„ calcu- 
lated using the Bondi rate as described in Sect. |4] and m <C 0.01, it was found that 
<~ 1 /5 of L acc is converted into jet power (|01)- The lowest m among the models in 
Fig.[5]is ~ 10~ 4 , and then accretion is highly radiatively inefficient, but the largest 
ism = 0.04, [...]. 
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For a representative model [obtained with the high resolution simulations with 
radiative and mechanical feedback of Ciotti et al. (2010), similar to the testcase 
ETG (Sect. 2.3.3)], the duty cycle (the fraction of the time the AGN is in the "on" 
state) is of the order of 10~ 2 , for the past 5-7 Gyrs; outside the nuclear bursts, the 
flow behavior is similar to that described in Sect.l3l( ll25Tl ). with a major difference: 
the lower central gas density, due to the MBH heating (13). [...] The brightness 
profiles are then much less centrally peaked than those in Sect. 2.3.3. [...] 

In addition to this important (positive) effect on the surface brightness profile of 
the hot gas, how does AGN feedback modify the scenario outlined in the previous 
Sections? ETGs where L$n > L~ v , and then already outflow-dominated, will not 
be affected by further sources of heating. For the other ETGs, the answer depends 
on how much energy from accretion is transferred to the hot ISM: if this energy is 
>> L$n then the scenario above will be modified, while if it is <<L$n it will be pre- 
served. In general, it can be noted that the gas modeling based on realistic stellar and 
dark mass profiles, stellar mass loss and supernova rates and their secular evolution, 
without accretion feedback can already reproduce reasonably well the fundamental 
gas properties (e.g., trend of Lx.ism with Lg, wide variation in Lxjsm, average gas 
temperature), therefore such modeling must catch the bulk of the origin and evolu- 
tion of the hot gas in ETGs. Moreover, even in the context of feedback modulated 
gas flow evolution, the hot gas content at the present epoch, seem still sensitive to 
the structural galaxy parameters, in the same sense as described in Sects. l2~3l and[3l 



(Ciotti & Ostriker, this volume; 19810 . Finally, the modeling without feedback - if 
any - shows the need for gas accretion from outside or confinement (Sect. [3}; the 
nuclear energy input should then mostly readjust the internal gas structure, without 
causing major degassing at later epochs. The measure in which activity affects the 
gas content is yet to be established observationally; so far, exploiting Chandra reso- 
lution, it has just been shown that the nuclear X-ray luminosities of ETGs correlate 
only weakly with their gas luminosity ( [199(1 ). 



5 Gas Flows and Galactic Shape, Rotation, Stellar Profile 

In Sect. [3] it was shown how the gas content of an ETG is sensitive to changes 
regarding the stellar and dark mass components that are in fact allowed for by ob- 
servations (see, e.g., the scatter around the fundamental scaling laws of ETGs), and 
by modeling (see, e.g., how mode l ETGs lying on the Fundamental Plane can be 
built with different <% and j3, lll05ll and Sect. 13.3. It . This holds even at fixed L B , so 
to account for a significant part of the large observed Lxjsm variation. Below we 
consider the effects on the hot gas content produced by additional variations in the 
galactic structure that are observed and have not been considered above, such as the 
galactic shape, the amount of rotation in the stellar motions, and the central stellar 
profile. 
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Fig. 8 The total X-ray emission (from ROSAT PSPC observations, (U) versus the galactic blue 
luminosity Lb (from total B-magnitudes in the Hyperleda catalog and distances in (§31) for a 
sample of ETGs with measured optical flattening and stellar rotation ( (9lll ). Galaxies with high 
flattening (b/a < 0.6) are marked with a red circle (see Sect. l5.1l for more details). 



5.1 Galactic Shape and Rotation 

Soon after the first large sample of ETGs with known X-ray emission was built from 
Einstein observations, it was found that the hot gas retention capability is related to 
the intrinsic galactic shape: on average, at any fixed Lg, rounder systems show larger 



total X-ray emission Lx and Lx/Lq than flatter elliptical and SO systems [39]. The 
relationship defined by Lx/Lq is stronger than that defined by Lx- Moreover, galax- 
ies with axial ratio b/a close to unity span the full range of Lx, while flat systems 
all have LxSjlO 41 erg s _1 (see, e.g., Fig. [8}. A similar result holds for the "diski- 
ness" or "boxiness" property of ETGs, that measures the deviation of the isophotal 
shape from a pure elliptical one (Ji, 64]). This property is described by the a^/a 



parameter, in a way that disky (an, > 0) ETGs show isophotes distorted in the sense 
of a disk, and boxy (04 < 0) ETGs have isophotal distortions in the sense of a box. 
Disky systems are also generally flattened by rotation, while boxy ones have various 
degrees of velocity anisotropy (see also Ir90lo . Boxy ETGs cover the whole observed 
range of Lx, while disky ETGs are less X-ray luminous on average (fl HH); this 
result is not produced only by disky galaxies having a lower average galactic lumi- 
nosity, with respect to boxy ETGs, since it holds even in the range of Lb where the 
two types overlap ( 19411 ). The relationship between Lx and cu/a was reconsidered, 
confirming the above trends, for the ROSAT PSPC sample ([138!]). 

There seems then to be a dependence of the hot gas content on the galactic shape, 
measured by either e = 1 — b/a or 04/0. Since flatter and disky systems also pos- 
sess, on average, higher rotation levels (J3l), the influence on the hot gas of both 
the shape of the potential well and of the stellar rotation was called into question. 
The gas in ETGs that are in part rotationally supported may have a lower "effec- 
tive" binding energy per unit luminosity compared to the gas in non-rotating ones 
(a lower effective L~ v /Lb ratio, in the notation of Sect. 12.31 see also Sect. 15.1.11 
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Fig. 9 The total X-ray emission versus the galactic luminosity Lg (both determined as in Fig.HJ; 
galaxies with v 10t /o" c > 0.5 are marked with a red symbol (see Sect. l5.1l for more details). 



below), and then rotating ETGs may be more prone to host outflowing regions. For 
this reason, the effects on Lx of the ellipticity e of the stellar distribution and of 
stellar rotation were studied for a sample of 52 ETGs with known Lx, maximum 
rotational velocity of the stars v rot , and central stellar velocity dispersion a c ( Il9lll ). 
[...]p The gas content can be high only for values of Vrot/Cc^O.4, while modest 
or low gas contents, as log[Lx(erg s _1 )/Lb(Lb i 0)]<3O.2, are independent of the 
degree of rotational support. Recently, for the ETGs of the SAURON sample, the 
relatio nship between soft X-ray emission and rotational properties was investigated 
again ( 111 1 ill ), confirming that slowly rotating galaxies can exhibit much larger lu- 
minosities than fast-rotating ones. As for the axial ratio and the isophotal parameter 
«4, the trend of Lx/L% with v Iot /(J c is not produced by the correlation between Lx 
and Lb : ETGs with high v ro t/Oc cover substantially the same large range in Lb as 
the whole sample (Fig. [9]). 

In conclusion, rotation seems to have an effect similar to that of shape, and Lx 
and Lx/Lb show a similar trend with respect to axial ratio, diskiness, and rotation: 
their variation is large for round, boxy and slowly-rotating systems, while it keeps 
below a threshold for flatter, disky and high-rotation systems. From observations it 
remains then undecided which one between axial ratio, diskiness, and rotation is 
responsible for the trend; more insights is given by the theoretical and numerical 
analysis discussed below. 



5.1.1 A Theoretical and Numerical Investigation 



The impact of stellar rotation and galactic shape on the hot gas content was also 
addressed with theoretical and numerical studies f ll3l l22l 13211 ). In principle, the 
lower gas content of flatter systems could be due to the mass distribution itself, or 
to a higher rotational level that decreases the effective potential. An analytic inves- 
tigation showed that flatter systems are less able to retain hot gaseous halos than 
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rounder ones of the same L^, due to the effect of the shape more than that of a larger 



rotational level ( 12210 . The investigation reconsidered the global estimate of the en- 
ergy budget of the gas introduced in Sect. 12.31 generalizing it for flows in flat and 
rotating galaxy models. The classical scalar virial theorem for a stellar distribution 
interacting with a dark matter potential can be written as 2T + 17 = \W\, where 71 
and T are the kinetic energies associated respectively with the stellar randorrQ and 
ordered motions, and \W\ is the potential energy of the stellar component plus the 
virial interaction energy of the stars with the dark halo. For a fixed total mass and 
mass distribution (i.e., a fixed \W\), the amount of rotational streaming energy T 
can formally vary from zero to a maximum that depends on the galaxy structure 



( 12210 : in the notation of Sect. 12.31 the power L rul related to rotational streaming 
is L m = M*T/M*, while that related to random motions is L a = 0.5M*7T/M*. 
How does L ro t enter the energy budget of the gas, for example in Eq. [9] for a fixed 
|W[ = 2M* (Lrot +L a ) /M* ? In an extreme case, the whole effect of the ordered 
motion is to produce a change in the effective potential experienced by the gas, 
if for example the galactic corona is rotating with the same rotation velocity as the 
stars; in this case L rot is to be subtracted fromL" v . In the opposite extreme case, all 
the kinetic energy of the gas, from random and from ordered motions, is eventually 
thermalized; then L rot is to be added to L$n and L a , in the denominator of Eq.[9] The 
real behavior, lying between the two extreme cases, can be parameterized re-writing 
Eq.|9]as 

L g r av -0CL w t 



L SN +L a + (l -d)Lr C 



with 0<a<l.lfa = 0, the thermalization of L rot is complete, and since the kinetic 
energy of stellar motions (L a +L rot ) will be lowei0than L$n, then Eq. [TTIcoincides 
with Eq.|9] If instead a = 1, there is no thermalization of L rot , the decrease of L~ rav 
is maximum, and the effect of rotation is maximum. However, it is found that the 
role of rotation remains minor, because it can change Eq.QjJby on ly a few P er cent: 
the variation of Eq. QT| between the null and the maximum L rot allowed by realistic 
galaxy models, is small, even for a = 1 (< 10%; |22To . Instead, variations of more 
significant amount that can make the gas significantly less bound (a decrease in L~ rav 
of ~ 20%) can be produced by a change in the galaxy structure, as a reasonable 
flattening of a round system at fixed Lg . Therefore, SOs and non-spherical ellipticals 
are less able to retain hot gaseous haloes than are rounder systems of the same Lb, 
and more likely to host outflowing regions. 

The results of the purely energetical approach above were tested with numerical 
studies of gas flows. Two-dimensional simulations for oblate ETGs, with different 
amounts of ordered and disordered kinetic energies, were carried out for gas in the 



inflow state ( 111310 . In this investigation Lx is reduced in rotating models, because 
the gas cools on a disk before entering the galactic core region, and then L+ v (Eq. O 
is reduced; since rotation increases on average with flatness, rotation would be the 



4 In the notation used here and in Sect. 12.31 77 is twice the energy due to random motions. 

5 For a totally velocity dispersion supported galaxy, L a accounts for the whole energy input to the 
gas from the stellar motions, that is significantly lower than Lsn (Sect. [531 below Eq.|5j. 
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underlying cause of the X-ray underluminosity of fiat objects. However, the massive, 
rotationally supported, and extended cold disk that forms in the equatorial plane, due 
to mass and angular momentum conservation, and comparable in size to the effective 
radius, is not observed; also, the resulting X-ray images should be considerably 
flattened towards the equatorial plane out to an optical effective radius or beyond, 
a phenomenon that is small or absent ([54]). Other authors ( 13211 ') performed two- 
dimensional numerical simulations of gas flows for flat systems, but allowing for the 
gas to be outflowing. The flows then developed a partial wind in flat ETGs that, if 
spherical, would be in inflow. In this way, the models accumulate negligible amounts 
of cold gas on a central disk. Rotation could also decrease the X-ray emission (of 
a factor of two or less), because it favoured the wind. In this scenario, then, flat 
models, rotating or not, can be significantly less X-ray luminous than spherical ones 
of the same Lb. because they are in partial wind when the spherical ones are in 
inflow; rotation has an additional but less important effect. 



5.2 The Central Stellar Profile 

[...] Interestingly, the radio luminosity Lr shows the same behavior as the total soft 
X-ray emission with respect to the inner stellar light profile: cusp ETGs are confined 
below a threshold in Lr, while core ones span a large range ofLje(J3[IS|9l).Core 
systems can then reach the highest Lr and possess a conspicuous radio activity cycle, 
while in cusp galaxies the radio emission keeps smaller, likely because of a rapid jet 
failure due to the lack of a dense confining medium, or a smaller duty cycle ( 19911 
and references therein). 
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